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2S HYPERFINE SPLITTING OF MUONIC HYDROGEN 

A.P.Martynenko* 

Samara State University, Theoretical Physics Department, 
^ ■ Pavlov Street 1, 44^011 Samara, Russia 

^ . Corrections of orders a , a are calculated in the hyperfine splitting of the 25 state 

CN i in the muonic hydrogen. The nuclear structure effects are taken into account in the 

one- and two-loop Feynman amplitudes by means of the proton electromagnetic form 
' factors. Total numerical value of the 25 state hyperfine splitting 22.8148 meV in 



X 



the fj,p can be considered as reliable estimation for the corresponding experiment 
with the accuracy 10~^. The value of the Sternheim's hyperfine splitting interval 
[8A£;^^"^(25) - AE^^^{1S)] is obtained with the accuracy 10"^ 



> 

O ■ PACS numbers: 31.30.Jv, 12.20.Ds, 32.10.Fn 

O 

^ I. INTRODUCTION 

The investigation of the energy levels of hydrogenic atoms with high accuracy is impor- 
tant for the check of the Standard Model and extraction the values of fundamental physical 
^ . constants (the fine structure constant, the electron and muon masses, the proton charge 
radius etc.) The Zemach radius takes substantial role among many fundamental param- 
eters determining the structure of the energy levels of simple atomic systems It is just 
this parameter which is responsible for the hyperfine splitting (HPS) interval in electronic 
and muonic hydrogen. The ground state hydrogen hyperfine splitting measurement was 
made many years ago with very high accuracy f3]: 

Ai/£J^(ep) = 1 420 405.751 766 7(9) kHz. (1) 

The theoretical expression for the hyperfine splitting in the hydrogen 

= E^(i + + 6''^ + 8^°' + 5^^^) , = 8 4 /^P^i^2 (2) 
^ ^3 (mi + m2j 

includes several corrections to the Fermi energy E^: S^^^ represents the QED contribution, 
the corrections 5'**'' and 5^°' are the proton structure and polarizability contributions, S^^^ is 
the contribution of hadronic vacuum polarization (HVP). fip is the proton magnetic moment 
in nuclear magnetons, mi is the lepton mass, m2 is the proton mass. The expression (2) is 
valid both for the muonic and electronic hydrogen but the exact value of the corrections is 
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TABLE I: The values of the Zemach radii. 



The Zemach 
radius 


Hyperfine splitting 
in electronic hydrogen [21] 


Hyperfine splitting 
in electronic hydrogen \22\ 


Elastic e-p 
scattering [23] 


Elastic e-p 
scattering [24] 


Rp 


1.047(19) fm 


1.040(16) fm 


1.086(12) fm 


1.013(16) fm 



essentially different for these atoms. The proton structure and polarizability effects lead to 
main theoretical uncertainty in the expression (2). The investigation of the nuclear structure 
and polarizability corrections in the h ype rfine splitting of the energy levels was done in many 
papers JM B H, lu IM 1^ E3, U, EM HE E3, EB, (see complete discussion in the review 
article |l7|). It is valid to say that during several tens of years the theory fall behind the 
experiment in solving this task. Existing difference between the theory and experiment 
without accounting the proton polarizability contribution can be expressed as follows [l7| : 



P) 



P) 



E^{e p) 



-4.5(1.1) X 10" 



(3) 



There exists at least two possibilities to improve theoretical expression (2). First possibility 
is connected with new more exact measurement of the electric G e and magnetic Gm proton 
form factors. The integral over the product of these form factors gives main part of the 
proton structure correction 6^^^ (the Zemach correction [3]): 



E 



p2fia f dp 



Ge{-p')Gm{-p' 



Hp 



E^{-2fia)Rp, W = afi, 



(4) 



where /i is the reduced mass of two particles, Rp is the Zemach radius. In the coordinate 
representation the Zemach correction (4) is determined by the contraction of the charge 
pE and magnetic moment pm distributions. The Zemach radius represents the integral 
characteristic of the proton structure effects in the hyperfine splitting of the energy levels. 
It may be considered as new fundamental proton parameter in the hydrogen atom. 

Another possibility consists in the using the muonic hydrogen. The measurement of 
the hyperfine splitting of the IS* and 25* energy levels with the accuracy 10~^ in the muonic 
hydrogen [18, 20., 21] allows to obtain the value of the Zemach radius with the accuracy 10~^. 
Then it can be used for obtaining new theoretical value for the hyperfine splitting in the 
electronic hydrogen and new theoretical restriction on the value of the proton polarizability 
contribution. An attempts to obtain the value of the Zemach radius were made recently 
using just as the experimental data on elastic electron-proton scatterin g so the hyperfine 



splitting in the hydrogen atom (see the Table 1 and the references 
difference between corresponding results reaches 7%. 
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The investigation of different contributions to the energy levels of the muonic atoms was 



done many years ago in Ref. jisl . 26l 27 1. So, at present there is need for new more complete 
analysis of all possible corrections in the HFS of the pp with the declared accuracy 10~^. The 
corrections of order to the hyperfine splitting of the 25" state and the Lamb shift 2P — 25* 
in the pp were calculated in Refs. 27, 2^|. In this study we continue the investigation of 
the different contributions of orders and to the muonic hydrogen HFS of the S states 
which was begun in Ref. |29| . The aim of the work consists in obtaining the numerical 
value of the 25* HFS with the accuracy 10~^ and the Sternheim's hyperfine splitting interval 
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FIG. 1: Effects of the one- and two-loop vacuum polarization in the one-photon interaction. 



in the muonic hydrogen with the accuracy 10 
serve as rehable guide for corresponding experiment and the measurement 2P 
shift |31|. Basic problems of the HFS measurement in the muonic hydrogen were discussed 



^ which can 
25* Lamb 



in Refs.|20|,|21 



II. EFFECTS OF VACUUM POLARIZATION IN THE ONE-PHOTON 

INTERACTION 

Our calculation of different energy levels of the hydrogen-like atoms are carried out on 
the basis of the quasipotential a ppr oach where the two-particle bound state is described by 
the Schroedinger-type equation [3^: 



where 



El - ml = El -ml 



''2 5 

/iR = E1E2/M is the relativistic reduced mass, M = Ei + E2 is the bound state mass. 
The quasipotential of the equation (5) is constructed in the quantum electrodynamics by 
the perturbative series using projected on positive states the two-particle off mass shell 
scattering amplitude T at zero relative energies of the particles: 

y = + \/(2) ^ y(3) ^ T = +T(2) +2-(3) ^ (6) 

We take the ordinary Coulomb potential as initial approximation for the quasipotential 
F(p, q, M): V{v. q, M) = V^{v - q) + A\/(p, q, M). 

The increase of the lepton mass when we change the electronic hydrogen to the muonic 
hydrogen leads to the decrease of the Bohr radius in the /ip. As a result the electron Compton 
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wave length and the Bohr radius are of the same order: 



n 



lie^ rrieC 



0.737384 (8) 



(me is the electron mass, is the reduced mass in the atom /ip). An important consequence 
of last relation is the increase the role of the electron vacuum polarization effects in the 
energy spectrum of the /ip ^33]. The effects of the vacuum polarization in the one-photon 
interaction are shown in Fig.l. 

To obtain the contribution of the diagram (a) Fig.l (the electron vacuum polarization) 
to the interaction operator there is need to make the following substitution in the photon 
propagator [33]: 

F 7r7o P(l-t;2) -4m2' ^ > 

At (— fc^) = ~ fj,l(Zay ~ ml{ZaY (electronic hydrogen, /ie is the reduced mass in 
hydrogen atom) we obtain — a/lSvrmg omitting the first term in the denominator of right 
part of Eq.(9). But when ~ fi'^{Za)'^ ~ m\(Za)^ (muonic hydrogen, m\ is the muon 
mass) than /ia and mg are of the same order and it is impossible to use expansion over a 
in the denominator of Eq.(9). To construct the hyperfine part of the quasipotential in this 
case (the muonic hydrogen) in the one-photon interaction we must use exact expression (9). 
We take into account that the appearance of the electron mass me in the denominator of 
the amplitude leads effectively to the decrease the order of the correction. It is well known 
that the hyperfine splitting quasipotential has the form |3^ : 

<^^(k) = ^i±^-l[(^,.,)k^ - (.ik)(.,k)]. (10) 

' m\m2 4 k^ 

For the S'-states 

SvrZ^^ (11) 
^ 3mim2 4 ' 

K=l. 792847337(29) is the proton anomalous magnetic moment. Averaging the potential 
(11) over the Coulomb wave functions we obtain main contribution of order (Za)^ to the 
HFS of the 25'-state in the system fip (the Fermi energy): 

1 ij^ 

E^{2S) = -{ZaY-^il + k) = 22.8054 meV. (12) 

The modification of the Coulomb potential due to the vacuum polarization (VP) is de- 
termined by means of Eq.(9) in the momentum representation as follows j3^: 

In the coordinate representation we obtain: 



VV^pW = ^J; d C' ( -— e--«^ ) . (14) 
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The contribution of the electron vacuum polarization to the hyperfine splitting part of the 
I7 quasipotential for the 5'-states can be derived in a similar way in the momentum and 
coordinate representations: 



3mim2 4 



a v/^^^(2^^ + 1) 
TT Ji 3^4(k2 + 4m2^2) 



3mim2 



TT Ji 



3^4 



TTOlr) e ^ 

r 



(15) 



(16) 



Using Eq.(16) we can obtain the electron vacuum polarization correction of order to 
the HFS in the np. Taking the wave function of the 25'-state 



-Wr/2 



1 



2V27r 

we represent this correction in the form: 



Wr\ 



W = fiZa, 



(17) 



AE- 



i7,yp 



HFS 



3mim2 IT 3W^ Jm^/w C 



1 



mi / 



(18) 



X 



1 



-r(25+l)/2C 



rdr 



0.0481 meV. 



The contribution of the muon vacuum polarization (MVP) can be found by means (16) after 
the substitution rrie — > mi. This correction is of order due to the reason mentioned above. 
Numerical value is equal 



3 fi 



AEr,l,yp = E'i2S)-^Za^ 

8 mi 



0.0004 meV. 



(19) 



In this expression and below E^{2S) is the Fermi energy of the 25* state. The diagrams of 
the two-loop electron vacuum polarization shown in Fig.l (b,c,d) give the contributions of 
the same order a^. The interaction operator corresponding to the loop after loop amplitude 
can be obtained using the relation (9). In the coordinate representation 



V"^^ (r 



8nZa(l + k) icricr2) 



3mim2 
°° Vr?2 - 1(27^2 + 1) 



3r/4 



5{v) 



m 



v^^{2e + 1; 



3^ 



X 



(20) 



Txr{rf — ^2) 



-2merir 



i -2m.eir 



and the contribution to the energy spectrum 

^-^17, VP-VP " 



0.0002 meV. 



(21) 



To calculate the contributions of the diagrams b, c in Fig.l which are determined by the 
polarization operator of the second order it is necessary to make the substitution in the 
photon propagator js^: 



1 



a 



f{v) 



TT/ Jo Ami + ^^(1 ~ 



dv 



a 



dv 



TT/ 3 7o 4m2 + A;2(i — t)2 



X 



(22) 



6 



Li2 + 2Lt2 ( — — I + - In z In — In ^ Inf 



l+V 



1 — V 



1 + vJ 2 1-v 



3 

-11(3 — f ^) In — ^ 211(3 — v'^) In 



1 -t; 



+ 



To find numerical value of this correction we write the quasipotential in the coordinate 
space: 

,2 



AyHFS ^ 
^^17, 2-loop VP) 



2 rl 



f{v)dv 



3mim2 \7r/ Jo (1— f 



6{r) 



7rr(l — v^) 



The potential (23) gives the contribution to the 25" HFS in the muonic hydrogen 



A^f'izoop yp = 0.0001 meV. 



(23) 



(24) 



We calculate all contribution numerically and the results are presented with the accuracy 
0.0001 meV. 



III. SECOND ORDER OF THE PERTURBATION THEORY 



The corrections of the second order of the perturbative series in the energy spectrum are 
defined by the reduced Coulomb Green function (RCGF) [361 : 



G2{r,r') =J2~gUr,r')YMyLi^')- 

Lm 



(25) 



The radial wave function gniir, r') was obtained in Ref. [2fl] as an expansion over the Laguerre 
polynomials. For the 25' - state 



m(m — 2) 



where x = fiZar, are the Laguerre polynomials: 



+ (^- + x-+x-je-^L}(x)Li(x^ 

(26) 



n\ \ dx 



dr. ( 



e 



(27) 



Some terms of the quasipotential contain the 5(r) so we have to know the quantity G2(r, 0). 
The expression for the RCGF was found in this case in Ref . js^] on the basis of the Hoestler 
representation for the Coulomb Green function after the subtraction the pole term: 



G25(r,0) 



Zan e 



2 „-x/2 



4x(x - 2)(lnx + C) + - 13x^ + 6x + 4 



(28) 



An 2x 

where C = 0.5772... is the Euler constant. The main contribution of order in the second 
order of the perturbation theory can be written in general form: 

< r2\vfpm X rn\^v,^^'\r2 > 

n=l,n^2 



A rpHFS _ 
^-^1 SOFT — 



^9 - E^. 



(29) 
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where /SV^^^ ~ 5{r). Using the relations (14), (28) we can present Eq.(29) as follows: 



a 



^Ers'oPT = E'{2S){l + a,)- I di 



(30) 



c/xe"^^^+^) 4x(x - 2)(lnx + C) + - 13x^ + 6x + 4 



0.0746 meV. 



The contribution of order in the second order of the perturbative series which is de- 
termined by the vacuum polarization can be derived from Eq.(29) changing ^V-^^^ — *• 
/SV^^yp. Using exact expressions for the wave function ip2{j) (17) and the RCGF (28) we 
write this correction 



^E^soPT = E^ (25') 



a 



97r2 

~ Vr/2 - 1(1 + 27^2) 



di]- 



(31) 



Fi(n = ^ -12 + 23ai - Sal - 4a? + 4af + 4ai(3 - 4ai + 2al) Inoi 



ai 



1 + 



2me 



2mii^ 1 

1^2 ^5^5 



(32) 

-1262(3 - 462 + 1h\) + 6a^(-4 + II62 - 862 + 262) - 4a^(-12 + 3962- 



-3462 + 126^) + 3a2(-12 + 5I62 - 5262 ^ 226^) , 02 = ai, 62 = ( 1 + -^r^ 



2mp 



2c2 



16mg^ 



as 



4aa 



r4m2 



as 



-21n(l-a3) 



+ 



+ 



1 - as (1 - as) 



-21n(l-a3) 



+ 



19as 



+ 



6ai 



ai(l + as) 



21n(l-a3)]}, 



16m2^?7 L 1 — as (1 — as)2 
1 



(l-«3)^ 

Numerical value of this contribution is equal 



as 



1 + 



w 

2m„l 



(1 



W 



^E^soPT = 0.0003 meV. 



(33) 



The second order of the perturbative series gives also other relativistic corrections of order 
(Za)^ including recoil effects which were studied in Ref. jl^jli^li^]. Corresponding numerical 
data are in the Table 1. 



IV. PROTON STRUCTURE AND VACUUM POLARIZATION EFFECTS 



The proton structure corrections in the system fip are relatively large in the comparison 
with the electronic hydrogen. In the HFS of the muonic hydrogen these corrections are 
defined in the leading order by the one-loop diagrams in Fig.2. 
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a b 

FIG. 2: Proton structure corrections of order (Za)^. Bold circle in the diagram represents the 
proton vertex operator. 

To construct the quasipotential corresponding to these diagrams we write the proton 
tensor: 



. 2m.2 



P2 - + "'2 



(p2 — k)'^ — ml + iO 



2m2 



(34) 

where p2, q2 are four momenta of the proton in initial and final states. The construction 
of the potential can be essentially simplified using the projection operators for the system 
muon-proton on the states with definite spin: 



(1+7°) 

2V2 



75, ^CSi) = [u{p2)v{pi)]s=i 



(1 + 7°) , 

2V2 



(35) 



where e'^ is the polarization vector of the state with the spin 1. Neglecting relative motion 
momenta of the particles in the initial and final states we obtain 



87r(l + fi;) '"7 7r2(A;2)2 

+16k'^klFl + l29>k%FiF2 + Uk^F^ + ^2k^Fl + Uk^FiF2 
Transforming the integration in Eq.(36) to the Euclidean space 

j d^k = 4:71 j k^dk j sin^ k^ = kcoscp, 



2 .2^^Fi-Mk%Fi+ (36) 

777,2 "^2 



{k^ - Amlkl){k^ - AmlklY 



(37) 



we make analytical integration over the angle and present the correction (36) as the 
one-dimensional integral over the variable k: 



AE- = -E-(2S)^.,„/;fvW. 



(38) 
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abed 

FIG. 3: Vacuum polarization and proton structure corrections of order a{Za)^. Dashed line in the 
diagram represents the Coulomb photon. 



V{k) 



+ 



mim2 (^^ _ m2)k{k + ^ Ami + k^) 



2U2 



+QAFiF2k'' + 16F|fc2 + 



32F^mik 
ml 



+ 



AF^k^ 

2 

mi 



AF^k^ 
ml 



+ 



(mi — m2)k{k + ^jAml + k"^) 



X 



128F^m''2 + 128FiF2m^ - IGF'^k'^ - QAFiF2k'^ - ASF^k"^ 

To cancel infrared divergence in Eq.(38) it is necessary to add the contribution of the iteration 
term of the quasipotential (10) in the HPS of the fip: 



iter,str 



3 mim2vr 



°° dk 



(39) 



where the angular brackets represent the averaging of the interaction operator over the 
Coulomb wave functions and the index HPS shows the hyperfine part in the iteration term of 
the quasipotential (8). The sum of the expressions (38) and (39) coincides with the result of 
Ref. [27[. The integration in Eqs.(38) and (39) was done by means of the parameterization of 
the proton electromagnetic form factors obtained from the analysis of elastic lepton-nucleon 



scattering Numerically the proton structure correction of order (Za)^ is equal 



AE 



HFS 

iter, str 



-0.1518 meV 



(40) 



Moreover, the effects of the proton structure must be taken into account carefully in the 
amplitudes of higher order over a shown in Pig. 3. 

The contributions of the diagrams (a) and (b) in Pig. 3 to the potential can be found as 
for the amplitudes in Pig.2. taking into account the transformation of one exchange photon 
propagator as in Eq.(9). Corresponding correction to the HPS of the energy level is equal 



^^stryP 



-E^{2S) 



Za 



2^ 



87r(l + /t) n Jo A;2(l 



1 v^{l-^)dv 



+ Ami 



dkVvp{k), 



(41) 



where the potential Vvp{k) differs from V{k) in the relation (38) only by the factor fc^. 
Despite of the finiteness of the integral (41) the amplitude terms of the quasipotential in Pig. 3 
(a), (b) must be completed by two iteration terms shown in Pig. 3 (c), (d). Pirst addendum 
< X X AVyp^ > of order a{Zay must be subtracted because the 27 amplitudes (a) 
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and (b) in Fig. 3 produce lower order contribution. Second term < Vyp x x V^^^ > 
which is also of order alZa)"^ has the structure similar to Eq.(29) of the second order of 
the perturbative series. The contributions of discussed iteration terms to the HFS of the fip 
coincide: 

^ESlyp^,,, = -2<V'xGfx AV«/^ >«''^= -2 < V{^p x x /XV""^ >'"'^= (42) 

meii Jo Jo K^[l — v^j + 1 
Numerical value of the proton structure and vacuum polarization effects of the 27 amplitudes 

A<|;L, + 2AESl'yp_,^,^. = -0.0026 meV. (43) 

Hadronic vacuum polarization contribution to the HFS in the fip was studied in Ref . f4l| . 
Here we present it in the different form using the expressions (38) and (41): 

.HFS _ TpFrnc\ «(^") r 



Dividing the integration range over s on the intervals where the cross section of the e~^e~ 
annihilation into hadrons (p(s) = cr^(e^e~ —>■ hadrons) /"isa^^) is known from the experi- 
ment i^l we can make numerical integration in Eq.(44). Numerical value coincides with the 
result obtained in Ref.|4l|: 

AEf^l = 0.0005 meV. (45) 



V. PROTON STRUCTURE EFFECTS, SELF ENERGY AND VERTEX 
CORRECTIONS OF ORDER a{Zaf 



There exists real number of important contributions of order which are presented in 
Fig. 4, 5. Radiative corrections of these amplitudes including recoil effects were studied earlier 
both in the Lamb shift and HFS of the hydrogen-like systems 17, 4^, 3l- Radiative photons 
were taken in the Fried- Yennie (FY) gauge jisl . 46l 47 1 where the mass shell amplitudes don't 
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contain infrared divergences. Infrared finiteness of the Feynman diagrams in this gauge 
gives the possibihty to make standard subtraction on the mass shell without introducing the 
photon mass. Let us consider radiative corrections which are determined by the self-energ 
insertions in the muon line. The renormalizable mass operator in the FY gauge is equal 



S^(p) = —{p — rnf' I dx 

TT Jo 



—3px 



m\x + (mf — p^) (1 — x) 



(46) 



Making the insertion (46) in the lepton tensor of the two-photon exchange diagrams and 
using the projection operators (35) we can construct the hyperfine splitting part of the 
quasipotential for the amplitudes in Fig. 4. In this case as before the vertex of the proton- 
photon interaction is determined by electric and magnetic form factors because the typical 
loop momenta are of order the muon mass. The contraction of the lepton and proton tensors 
over the Lorentz indices and the Dirac 7 matrix trace calculation were made in the system 
Form [4^. In the Euclidean space of the variable k we can present the correction to the 
HFS of the muonic hydrogen as follows: 



A W^FS 



VsE{k,(j),x) 



STT^ TT Jo 



1 roo rn 

xdx / kdk / sin^ (j)d(pVsE{k, (p, x) 
Jo Jo 



X' 



[k"^ + 4m2 cos^ 0) [{xm\ + xk'^Y + 4mf x^A;^ cos^ 1 



5m 



X 



(47) 
(48) 



rF|(x + 6^) cos^ (p ^k'xF;^ + 16mfF2 cos* 0(4Fix - F2X - 2F2x) + 



ml 



+ 16ml cos^ 4>{F^x + QFlx + AF1F2X + 8F1F2X + F|x + 2F|x) + 



Ak 



8k^ 



+32mixFi{Fi + F2) jF^x cos^ jF^x - IGk^F^x cos* 0+ 



rrio 



+16k^x cos^ (f){F^ + 4F1F2 + F|) + 32k^Fix{Fi + F2)|. 

After analytical integration over the angle we present the contribution (47) in the integral 
form which was used for numerical calculation: 



Qm\k^F2X 



+ 

2km\ 
ml 



4 

ml 



xdx 
Ak 



dk< 



-^ + 32Fi(Fi + F2) 
mi 



1 



hi{k, x) 



F2(2Fi + F2)x 



m2\xmi + xk"^) m2 
kF^ 



k 



(49) 



ml 



■{xm^ + xk ) 



\h2{k,x) hi{k,x) 
e ( 1 



+ 



lilik.x) ml{xmi + xk'^) \h2{k,x) 
/ii(/c, x) = kJ Amlx'^k'^ + (xmf + xk'^y + {xm^ + xk'^)\l Ami + k'^i 



hi{k, x) 



h2{k, x) = \J Am\x'^k'^ + (xmf + xk'^Y + {xm\ + xk"^). 
Numerical value is equal 



A T^HFS 
^-^27, 



0.0010 meV. 



(50) 
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abed 

FIG. 5: Proton structure and muon vertex effects of order a{Zctf' . 

Let us consider calculation of the vertex corrections. The renormalizable expression of 
the one-particle vertex operator in the FY gauge was obtained in Ref. j49| (pi = m^): 



a 



»1 rl 

dx dz 

4:71 Jo Jo 



F(2) 



A 



A2 



(51) 



where A = mix + 2pk{l — x)z — k'^z{l — xz), the functions Fj^), Fj^) were determined in 
Ref. li^ . The lepton tensor can be divided into two parts: 



'^{Pi)P!)^K~Pi — k + rni)'^^v{qi){k'^ — 2k'^mi)[m\x — k'^z{l — xz) + 2mik'^x'^] 
{k^ — AkQrn{)[{m'{x — k'^z{l — xz))^ — 4mf /EgX^-z^] 



(52) 



MW(2) 



v{pi)F^'^\—pi — k + mi)7^f (gi)(A;2 — 2k^mi)[mlx — k'^z{l — xz) + 2mik'^x'^]'^ 
(/c^ — 4/comf)[(mfa; — k'^z{l — xz))"^ — Arniklx^ z'^Y 



(53) 

Keeping for the simplicity the main contribution over mi/ 1112 we write this type vertex 
corrections as follows: 



A WHFS 



-E^{2S) 



a;\2 8mim2 '"^ '"^ 



TT 



dx I dz TTsin^ 

1 + KlTT Jo JO JO 



kdk X 



(54) 



X 



Vi{x,k,^)[Fi{Fi + F2)-{l + K)] 



2^21 



(/c^ + Ami cos^ (f)){k'^ + Am2 cos^ 0) [[mfa; + k'^z{l — zx)]'^ + Amlk^ cos^ (px 
Vi{x, k, 0) = -2m^x2(l -x) + Pml{6x^z^ - Sx^z"^ - Sx'^z + 8xz - 3x) + (55) 
+k^{4:x^z^ - Gx'^z^ - 5x^z^ + 12xz^ - 2xz^ - 6z^ + 3z), 



^EZlrt 2 = -E'{2S) (^)'^^ J\{l-x)dx I 'dz / vrsin^ 



k^dk X (56) 



X ■ 



V2ix,k,(P)Fi{Fi + F2) 



{k"^ + 4mf cos^ (j)){k'^ + 4m2 cos^ cj)) [[mfx + A;22;(l — zx)Y + Am^k"^ cos^ ^x^Z' 



t2 ' 
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V2{x, k, (j)) = mlx^z{2z - 1) - k'^mlxz^{Axz^ - 2xz - 4z + 2) + 
+k^z^{2x^z^ - x^z^ - Axz^ + 2xz + 2z- 1). 



(57) 



The iteration contribution is equal 



iter, 27 vert 



(58) 



After analytical integration in Eqs.(54) and (56) over the angle and the subtraction (58) 
(one photon is the Coulomb-like and the other one contains the hyperfine part of the potential 
with the value of magnetic form factor at zero point) we have the expressions of the diagrams 
(a) and (b) in Fig. 5: 



^-^^27, vert 



-E^(25) 



a 



dx I dz 







dk< 



Fi{Fi + F2) 



8A;(1 + K)'m\m2x'^ z"^ 



—2m\x'^x + k^mfx x 



2_2, 



(59) 



X (6x2^2 



ixz 



3xz + 8^ - 3) + k^ziAx-^z-^ - Qx^z^ - hx^z^ + I2xz^ - 2xz - 62 - 3) 



X 



+ 



+ 



+ 



Fi(Fi + F2)x 
2(1 + K,)mim2kx^z^ 



h{a? - }PW - c^) a(a2-62)(a2_c2) ^(6^ - c2)(a2 - c2)_ 
m\x'^z{2z - 1) - 2k'^m\xz^{2xz^ - xz - 2z + I) + A;V(2x^/ - x'^z^ - Axz^ + 2xz + 2^-1) 

VTT¥ ^/TT^ 1 



2c3(62-c2)(a2-c2) C(62 -c2)(a2 -c2) c(62 - c2) (a2 - c2)2 

2 /c^ 2 2 [mix + k'^z{l - xz)]"^ 



fc2 



Ami 



2 ' 



Ami ' 



Am\k'^x'^z'^ 



Numerical value of the vertex correction (59) is equal 

^E^-y^vert = -0.0018 mcV 



(60) 



(61) 



Next vertex type diagram with one rounded photon and two exchanged photons is the 
diagram of the "jellyfish" type. Its contribution to the energy spectrum is of order a{Za)^. 
At small loop momenta this diagram gives the finite answer in the FY gauge. The lepton 
tensor relating to the diagrams (c) and (d) in Fig. 5 was obtained in Ref . |i44] : 



a 



Att Jo 



xdx 



z)dz ^ 



71=1 



(62) 



where A has the form as in Eq.(51). The tensor functions M^"-* are written explicitly 
in Ref. 3|- The character of further transformations of the amplitudes (c), (d) in Fig. 5 
to construct the HFS part of the potential is the same as for other amplitudes shown in 
Figs. 4, 5. Omitting the details of such transformations which were carried out by means of 
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analytical system Form j48| we write here three contributions to the HFS corresponding to 
the functions in the leading order over the ratio rrii/m^. 

^Ef,feUyf.sH = -^M!/^ ^^^^ _ ^^^^ _ 3^^^ ^^^^^^^^ ^ ^ (gg^ 

f"^ siv? (f)d(f) [mix + k'^z{l — xz)] 







/o (fc^ + 4m2 cos^ 0) [mf a; + /c^2;(l — xz)]^ + 4mf /c^ cos^ 0x^2;^ ' 
A^S,,.. = _ ^^^^^^^^ ^ ^ 

sin^ (j)d(j) [mix + A;^2;(l — xz)]'^[k^xz^{l — xz) + m1{x^z + 2x2; — x — 3^;)] 
+ 4mi cos2 0) {[mfx + k^z{l - xz)]'^ + Amjk^ cos^ (f)x^z^}^ ' 

^E^,feUyf^sH = ^^"^^^^^ £ xdx j\l - z)zHz f k' dkmlF,{F, + F,) X (65) 
2 r sin^0(i0 [mix + k'^ z{l — xzyf 

The integration over the angle can be done in Eqs.(63)-(65) analytically. Omitting inter- 
mediate expressions we can write final numerical result to the HFS of the fip: 

= E ^E^,%yfisk = 0.0005 meV. (66) 

n=l 

In the point-like proton approximation when the nucleus form factors entering the Feynman 
amplitudes in Fig. 4,5 are changed on their values at fc^ = (-Fi(O) = 1, -^2(0) = k) the 
contributions (63)-(65) will increase twofold. 



VI. CONCLUSION 

The calculation of different quantum electrodynamical effects, effects of the proton struc- 
ture and polarizability, the hadron vacuum polarization to the HFS of muonic hydrogen is 
performed in this work. The corrections of orders and are considered. Working with 
the vacuum polarization diagrams we take into account that the ratio fia/me is very close to 
1 and don't increase the order of corresponding contributions. Obtained numerical results 
are presented in the Table 1. We include here also QED corrections to the Fermi energy 
which are determined by the muon anomalous magnetic moment a^E^{2S) ^3 (experi- 
mental value of the muon anomalous magnetic moment a^^^ = 11659203(8) x 10~^° p(f\ is 
used), the Breit relativistic correction of order (Za)^ jH^I, the proton structure corrections 
of order (Z a)^ \n(Z a)"^ M, t he hadron vacuum polarization contribution 41] and the proton 
polarizability correction [16[ , the weak interaction contribution due to the Z boson exchange 

Let us point out some peculiarities of this investigation. 

l.The effects of the vacuum polarization play very important role in the case of the 
muonic hydrogen. They lead to essential modification of the spin-dependent part of the 
quasipotential of the one-photon interaction. 
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TABLE II: Corrections of orders a^, to the 25" state HFS in the muonic hydrogen. 



Contribution to HFS of fj,p 


Numerical value in meV 


Reference 


The Fermi energy E^(2S) 


22.8054 


[17], (12) 


Muon AMM correction a^E^{2S) of order a^,a^ 


0.0266 


\n] 


Relativistic correction ^[Z of' [2S) of order 


0.0026 


[51] 


Relativistic and radiative recoil corrections 
with the account proton AMM of order 


0.0018 


[38] 


One-loop electron vacuum polarization 
contribution of I7 interaction of orders a^, 


0.0482 


(18) 


One-loop muon vacuum polarization 
contribution of I7 interaction of order 


0.0004 


(19) 


Vacuum polarization corrections of orders , 
in the second order of perturbative series 


0.0746 


(30)+(33) 


Proton structure corrections of order 


-0.1518 


[27], (40) 


Proton structure corrections of order 


-0.0017 


[7] 


Electron vacuum polarization contribution-l- 
proton structure corrections of order 


-0.0026 


(43) 


Two-loop electron vacuum polarization 
contribution of I7 interaction of order 


0.0003 


(21)+(24) 


Muon self energy -|- proton structure 
correction of order 


0.0010 


(50) 


Vertex corrections -|- proton structure 
corrections oi orcier a 


-U.UUio 


(olj 


"Jellyfish" diagram correction -|- 
proton structure corrections of order 


0.0005 


(66) 


HVP contribution of order 


0.0005 


(45) 


Proton polarizability contribution of order 


0.0105 


[16] 


Weak interaction contribution 


0.0003 


[52] 


Summary contribution 


22.8148 ± 0.0078 





2. The proton structure effects are taken into account consistently in the loop ampli- 
tudes by means of electromagnetic form factors. The point-like proton approximation gives 
essentially increased results (approximately twofold). 

3. The calculation of the muon self-energy and vertex corrections of order a{Za)^ is 
done on the basis of the expressions for the lepton factors in the amplitude terms of the 
quasipotential obtained by Eides, Grotch and Shelyuto. We supplemented these relations 
by the subtraction of the iteration terms of the potential. 

Total value of the 25* state HFS in the muonic hydrogen shown in the Table 1 can be 
considered as definite guide for the future experiment which is prepared Numerical 
values of the corrections were obtained with the accuracy 0.0001 meV. Theoretical error 
connected with the uncertainties of fundamental physical constants (fine structure constant, 
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the proton magnetic moment etc.) entering the Fermi energy compose the value near 10~^ 
meV. Other source of theoretical uncertainty is connected with the corrections of higher 
order. Its estimation can be found from the leading correction of the next order on a and 
mi/m2 in the form: a{Z aY \\i{Z aY {2S) / ti ^ 0.00003 meV (the value of fine structure 
constant is a = 137.03599976(50) Q). Further improvement of theoretical result presented 
in the Table 1 is connected first of all with the corrections on the proton structure and 
polarizability which give the theoretical error near 34 x 10~^ ppm. The most part of this error 
is determined by the proton structure corrections of order {ZaY (the Zemach correction). 
So, the measurement of the hyperfine splitting of the levels 15" and 25* in the muonic hydrogen 
with the accuracy 30 ppm will lead to more accurate value (with relative error 10~^) for 
the Zemach radius which than can be used for the improvement of theoretical result for 
the ground state hydrogen hyperfine structure and more reliable estimation of the proton 
polarizability effect. 

Another important quantity regarding to the hyperfine structure of the muonic hydrogen 
is the Sternheim's hyperfine splitting interval [3u\. It doesn't contain the proton structure 
and polarizability effects of the leading order and provides additional test of quantum elec- 
trodynamics for the hydrogen atom. Accounting numerical results of this work and Ref. |2j| 
we can write numerical value for this interval in the form: 

[8AE^^^(25) - AE^^^(IS)] = -0.120 meV. (67) 

This result is valid with a precision 10^^ due to the cancellation of the proton structure 
and polarizability corrections. The increase the number of the tasks due to excited states of 
simple atoms [53| and the inclusion new simple atoms where the hyperfine structure of the 
energy spectrum is studied will decrease the uncertainties in the determination of physical 
fundamental parameters and increase the accuracy for the check of the Standard Model in 
low energy physics. 
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